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1. Introduction 



It is expected that strongly interacting matter undergoes a transition in some temperature in- 
terval from hadron gas to deconfined state also called the quark gluon plasma (QGP) Creating 
and exploring deconfined medium in the laboratory is the goal of the large experimental heavy 
ion program at RHIC [§] and LHC [||]. Studying hot and dense strongly interacting matter is also 
subject of a large effort in lattice QCD (see Refs. [Q, for recent reviews). Early lattice QCD 
simulations at non-zero temperature were limited to large quark masses and had no control over 
the discretization errors [|[ [7J §]. During the past seven years calculations with the physical strange 
quark mass and physical or nearly physical light (u,d) quark masses have been performed using 
improved staggered fermion actions 10 , [TT , 12, |13], 14[ 15,1^, [l^, 18, 19] and for several quan- 
tities continuum extrapolated results have been obtained. There was also progress in lattice QCD 
calculations at non-zero temperature using other fermion formulations, namely Wilson fermions 
[20, 21, 22], Domain-Wall fermions [23] and overlap fermions [24]. The later two formulations 
preserve the chiral symmetry of continuum QCD. However, due to much larger computational costs 
of these formulations the corresponding results are far less extensive. 

To get reliable predictions from lattice QCD the lattice spacing a should be sufficiently small 
relative to the typical QCD scale, i.e. AgcD« <C 1. For staggered fermions, discretization errors 
go like ^((aAgco) 2 ) but discretization errors due to flavor symmetry breaking turn out to be quite 
large numerically, and dominate the cutoff dependence of thermodynamic quantities at low tem- 
peratures. To reduce these errors one has to use improved staggered fermion actions with so-called 
fat links [25]. At high temperature the dominant discretization errors come from the lattice distor- 
tions of the quark dispersion relation and go like (aT) 2 , and therefore could be very large. Thus 
it is important to use improved discretization schemes, which reduce or eliminate these discretiza- 
tion errors. Staggered fermion actions used in numerical calculations typically implement some 
version of fat links as well as improvement of quark dispersion relation and are referred to as p4, 
asqtad, HISQ and stout. Independently of specific improvement all lattice results eventually must 
be extrapolated to the continuum limit. 

In this contribution I am going to discuss lattice QCD calculations on the equation of state, 
study of deconfinement aspects of the QCD transition, including color screening and fluctuations of 
conserved charges, determination of the chiral transition temperature. I will mostly discuss lattice 
results obtained with staggered quark formulation, where appropriate results obtained using other 
actions will be highlighted. 



2. Equation of State 

The equation of state has been calculated using different improved staggered fermion actions 
p4, asqtad, stout and HISQ. In the lattice calculations of the equation of state and many other 
quantities the temperature is varied by varying the lattice spacing at fixed value of the temporal 
extent N t . The temperature T is given by the lattice spacing and the temporal extent, T = 1/ (N T a). 
Therefore taking the continuum limit corresponds to N t — > °° at the fixed physical volume. The 
calculation of thermodynamic observables proceeds through the calculation of the trace of the 
energy momentum tensor, s — 3p, also known as the trace anomaly or the interaction measure. 
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This is due to the fact that this quantity can be expressed in terms of expectation values of local 
gluonic and fermionic operators, (see e.g. Ref. [|5|]). Different thermodynamic observables can be 
obtained through integration of the trace anomaly 1 . The pressure can be written as 

P( T) m_r*- le _ 3p) . (2 „ 



r 4 r 4 Jt t' 5 

The lower integration limit To is chosen such that the pressure is exponentially small there. Fur- 
thermore, the entropy density can be written as s = (e + p)/T . Since the interaction measure is 
the basic thermodynamic observable in the lattice calculations it is worth to discuss its properties 
more in detail. In Fig. [I] (left panel) I show the results of the calculation with different improved 
actions. The calculation with p4 and asqtad actions have been performed on N T = 8 lattices and 



light quark masses to/ = m s /10, with m s being the physical strange quark mass [13, 15]. These 



light quark masses correspond to the pion masses slightly above 200 MeV in the continuum limit. 



For this value of N r the above deviation from the physical quark mass plays little role [14, 26]. The 
N T = 12 asqtad calculations have been performed for to/ = m s /20 p7[]. Calculations with HISQ 
action have been performed for N t = 4, 6, 8, 10 and 12 for to/ = m s /20 corresponding to the 
pion mass of 160 MeV in the continuum limit [27]. A subset of these results is shown in Fig. |]. 
Finally, calculation of the trace anomaly and the equation of state was performed with stout action 
using N t = 4, 6, 8, 10 and 12 and physical light quark masses [111]. Using the lattice data from 



N T = 6, 8, 10 a continuum estimate for different quantities was given J18Q. The interaction mea- 
sure shows a rapid rise in the transition region and after reaching a peak at temperatures of about 
200MeV decreases. At low temperatures lattice data obtained with HISQ action and stout action 
agree with each other. We also compare the lattice results with hadron resonance gas (HRG) model 
calculations which seems to agree well with HISQ and stout results for T < 150 MeV. Cutoff ef- 
fects (i.e. N t dependence) appears to be the strongest around the peak region. They decrease at 
high temperatures and at T > 400 MeV all lattice results agree with each other. At low temperature 
the cutoff effects related to flavor symmetry breaking are very large for p4 and asqtad actions. Due 
to these large cutoff effects the N T = 8 p4 and asqtad data are below the hadron resonance gas 
model. Taking into account the distortions of the hadron spectrum due to flavor symmetry break- 
ing in hadron resonance gas calculations leads to good agreement of HRG model with the lattice 



[]28[]. Since at high temperatures the effects of flavor symmetry breaking in the pressure and the 
interaction measure are small the reduction in e — 3p at low temperatures must be at least partly 
compensated by a larger value at intermediate temperatures. Thus, the large N T dependence of the 
peak height of e — 3p is related to the large flavor symmetry breaking effects for p4 and asqtad 
actions. 

In Fig. [I] I also show the entropy density divided by the corresponding ideal gas value and 



compare the lattice results with resummed perturbative calculation [ |29| , p0| ], as well as with the 
predictions from AdS/CFT correspondence for the strongly coupled regime [31]. The later is con- 
siderably below the lattice results. Note that the pressure, the energy density and the trace anomaly 
have also been recently discussed in the framework of resummed perturbative calculations which 



seem to agree with lattice data quite well at high temperatures [p2|] . The differences between the 



A somewhat different approach was used in Ref. [18| 
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Figure 1: The interaction measure (left) and the entropy density (right) as function of the temperature 
calculated with improved staggered fermion actions. The dashed line in the left panel shows the HRG result. 
The band in the right panel shows the resummed perturbative result, while the solid line is the expectations 
based on the strongly coupled limit. 



stout action and the p4 and asqtad actions for the trace anomaly translates into the differences in 
the pressure and the entropy density. In the high temperature region the entropy density obtained 
with stout action is 10% smaller than the entropy density obtained with p4 and asqtad actions. 

3. Taylor expansion of the pressure and fluctuations of conserved charges 

Due to the infamous sign problem lattice QCD Monte-Carlo simulations are not possible at 
non-zero quark chemical potentials. The pressure and other quantities at non-zero chemical po- 
tentials, however, can be evaluated using Taylor expansion. The Taylor expansion can be set up 
in terms of the quark chemical potentials \i u , and or in terms of the chemical potentials 
corresponding to baryon number B, electric charge Q and strangeness S of hadrons 

P 1 , .. .. N V- 1 



T 4 y T 3 



v uds = : r 1 : n ii 

d(n u /Tyd(n<i/Tyd(n s /T) k 1 " ; 



d i+ J +k p/T 4 



P_ 

j, 4 



d(n B /Tyd(nQ/T)jd(n s /Tr K ■ } 

Using Taylor expansion method the equation of state has been calculated for small chemical po- 



tential in the continuum limit Q33p . Earlier results at non-zero lattice spacing have been reported in 
Refs. [34, 35]. While Taylor expansion can be used to study the physics at non-zero baryon density 
the expansion coefficients are interesting on their own right as they are related to the fluctuations 
and correlations of conserved charges. As will become clear later fluctuations and correlations of 
conserved charges are good probe of deconfinement. Fluctuations of conserved charges are also 



useful for determining the freeze-out conditions in heavy ion experiments [36]. 
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Figure 2: The fluctuations of baryon number (left), electric charge (middle) and strangeness (right) as 
function of the temperature calculated with HISQ and stout action in the continuum limit and normalized by 
the corresponding ideal quark gas values % 2 ■ The solid black curves correspond to HRG predictions. 



The diagonal expansion coefficients are related to second and higher order fluctuations of 
conserved charges 

ll = ^(N 2 x) 

X? = ^(("x)-HN 2 x) 2 ) (3-3) 
while off-diagonal expansion coefficients are related to correlations among conserved charges, e.g. 

Xu=^(NxN Y ). (3.4) 

Second order fluctuations have been studied with improved staggered actions in Refs. [15, 17, p7| ]. 
Recently continuum results have been obtained for second order fluctuations of baryon number, 



electric charge and strangeness using stout and HISQ actions [ |38| , |39| ] which are shown in Fig. 
^| The lattice results are also compared with HRG model. As one can see from the figure at low 
temperatures the lattice data are described well by HRG model indicating that the dominant degrees 
of freedom in that temperature range are hadronic. Deconfiniment is seen as a rapid increase of the 
fluctuations for T > 150 MeV, which eventually reach values that are close to the expectations of 
weakly interacting quark gas. In other words, fluctuations indicate quark degrees of freedom at high 
temperatures. The lattice results obtained with HISQ and stout action agree well in the continuum 
limit, except for temperatures around 200 MeV where some discrepancies are observed in baryon 
and electric charge fluctuations. Another way to study deconfinement is to consider correlations 
of conserved charges. These correlations are very different for hadron gas and quark gas. As an 
example let us examine the strangeness-baryon number correlation. It is convenient to normalize 
this quantity as follows Cbs = —^xfi/Xs- At high temperatures where strangeness is carried by s- 
quarks this quantity should be close to one. At low temperatures on the other hand strange baryons 
are responsible for strangeness-baryon correlations. In Fig. ||(left) I show continuum results for 



Cbs obtained with stout [38] and HISQ [39] actions. At low temperatures the lattice results are 
described by HRG model, while at high temperatures are close to one as expected for the quark 
gas. 

Fourth order fluctuations for baryon number, electric charges and strangeness have been stud- 



ied using p4 [37] and HISQ Q40p actions. Cutoff effects are quite large for p4 action in the low 
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Figure 3: Baryon number-strangeness correlation Cbs (left) and fourth order fluctuations of baryon number 
and electric charge for HISQ action normalized by the corresponding ideal quark gas value as function 
of T/T c (right). The lines correspond to the prediction of the hadron resonance gas model. Open (filled) 
diamonds in the right panel correspond to N T = 4 p4 results for %2 (X4 )■ 



temperature and the transition regions. These large cutoff effects also result in much larger value 
of the transition temperature than obtained with stout and HISQ actions [10]. Calculations of the 
fourth and higher order fluctuations is quite demanding computationally and for this reason no con- 
tinuum results have been obtained yet. The lattice results obtained with HISQ action for N T = 6 
and 8 are shown in Fig. || as function of the temperature in units of the chiral transition tempera- 
ture T c = 154 MeV (see below) and also compared with earlier results obtained with p4 action and 
N r = 4. For the p4 action we use the value of the transition temperature T c = 204 MeV determined 
in Ref. Jl0| ] for N z = 4. To reduce the cutoff effects the lattice spacing in the HISQ calculations 
was set by the kaon decay constant fg. At high temperatures the fourth order fluctuations are close 
to the values corresponding to non-interacting quark gas. At low temperatures the fourth order 
baryon number fluctuations are reasonably well described by hadron resonance gas. This is not the 
case, however, for the electric charge. One possible reason for the disagreement between the lattice 
and HRG model in this case could be the large cutoff effects in the pion sector. Electric charge 
fluctuations are very sensitive to the pion sector which is largely distorted on the lattice even if 
HISQ or stout action is used. These distortions correspond to effectively larger pion mass that 
would explain why the lattice data are below the HRG expectations. Fourth order fluctuations have 
a maximum in the transition region. Interestingly enough the position of the maximum in T /T c is 
roughly the same for the HISQ and p4 actions. The height of the peak, however, is significantly 
larger for the %f when p4 action is used. 

Second and fourth order light and strange quark number fluctuations at high temperatures have 
been studied on the lattice in Refs. [38, |d] and compared with resummed perturbative calculations 
[4-2, 43]. The resummed perturbative calculations seem to describe the lattice data quite well at 
T > 300 MeV. 



4. Deconfinement : color screening 

One of the most prominent feature of the quark gluon plasma is the presence of chromoelectric 
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(Debye) screening. The easiest way to study chromoelectric screening is to calculate the Polyakov 
loop. The Polyakov loop is an order parameter for the deconflnement transition in pure gauge 
theory, which is governed by Z(N) symmetry. For QCD this symmetry is explicitly broken by 
dynamical quarks. There is no obvious reason for the Polyakov loop to be sensitive to the singular 



behavior close to the chiral limit although speculations along these lines have been made [|45|]. The 
Polyakov loop is related to the screening properties of the medium and thus to deconflnement. After 
proper renormalization, the square of the Polyakov loop characterizes the long distance behavior 
of the static quark anti-quark free energy; it gives the excess in free energy needed to screen two 
well-separated color charges. The renormalized Polyakov loop, calculated on lattices with temporal 
extent N x , is obtained from the bare Polyakov 

L ren {T) = z(pf*L bare (p) = z(Pf* ( ^W(x) ) , W(x) = n U (x Q ,x), (4.1) 

where Uq = exp(jgaAo) denotes the temporal gauge link and z(/3) is the renormalization constant 



determined from the T = static potential Ql 1Q. Continuum results for the renormalized Polyakov 



loop have been obtained with stout [17] and HISQ actions [48]. These are shown in Fig. |] together 
with N r = 6 HISQ results [|l9|]. One can see a good agreement between the stout and HISQ results. 
I also compare the 2+1 flavor QCD results with the corresponding results in pure gauge theory p6[ 
47] as well as with the prediction of non-interacting gas of static-light(strange) hadrons p9| ]. 



We see that in the vicinity of the transition temperature the behavior of the renormalized Polyakov 
loop in QCD and in the pure gauge theory is quite different. The calculation of L ren based on non- 
interacting static-light hadron gas can explain the lattice data for T < 140 MeV. The renormalized 
Polyakov loop has also been calculated using lattice fermion formulations other than staggered, 
namely the Wilson formulation [^] and the overlap formulation [f£4]]. These formulations are 
considerably more expensive computationally than the staggered formulation and therefore the 
calculations have been performed at unphysical pion mass. The results obtained using Wilson 
action and overlap action for the Polyakov loop agree very well with the staggered fermion results 



at the same value of the pion masses fl2Q, 24]. 



Further insight on chromoelectric screening can be gained by studying the singlet free energy 



of static quark anti-quark pair (for reviews on this see Ref. (pQ, 51]), which is expressed in terms 



of correlation function of temporal Wilson lines in Coulomb gauge 

exp(-F 1 (r,r)/r) = l -Tv{W(r)W\0)). (4.2) 

Instead of using the Coulomb gauge the singlet free energy can be defined in gauge invariant 
manner by inserting a spatial gauge connection between the two Wilson lines. Using such definition 



the singlet free energy has been calculated in SU (2) gauge theory [|52|]. It has been found that the 
singlet free energy calculated this way is close to the result obtained in Coulomb gauge [p2|]. The 
singlet free energy turned out to be useful to study quarkonia binding at high temperatures in 
potential models (see e.g. Refs. [0]). It also appears naturally in the perturbative calculations of 



the Polyakov loop correlators at short distances [53] 



The singlet free energy has been recently calculated in 2+1 flavor QCD with HISQ action on 
24 3 x 6 and 16 3 x 4 lattices [54]. The numerical results are shown in Fig. [| At short distances the 
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Figure 4: Left: The renormalized Polyakov loop as function of the temperature in 2+1 flavor QCD and 
pure gauge theory. Right: The singlet free energy as function of the distance at different temperature calcu- 
lated with HISQ action [p4{], The solid line in the left panel corresponds to static-ligh(strange) hadron gas 
prediction for L re „ (see text). The dashed line in the right panel is the T = potential [n9h. 



singlet free energy is temperature independent and coincides with the zero temperature potential. 
In purely gluonic theory the free energy grows linearly with the separation between the static quark 
and anti-quark in the confined phase. In presence of dynamical quarks the free energy is saturated 
at some finite value at distances of about 1 fm due to string breaking (see e.g. Ref. [50]). This 
is also seen in Fig. ^. Above the deconfinement temperature the singlet free energy is exponen- 
tially screened at sufficiently large distances [46, 47] with the screening mass proportional to the 
temperature , i.e. 

F,(r,T)=F ca (T)-^^-ex V (-m D (T)r), m D ~T. (4.3) 

The lattice data for the singlet free energy are consistent wit these expectations for r > 0.S/T. 

Let me finally note that contrary to the electro magnetic plasma the static chormomagnetic 
fields are screened in QGR This is due to the fact that unlike photons gluons interact with each other 
(the stress tensor is non-linear in QCD). Magnetic screening is non-perturbative, i.e. it does not 
appear at any finite order of pertubation theory. In lattice calculations chromomagnetic screening 
is studied either in terms of spatial Wilson loops [^5j or in terms of spatial gluon propagators 



Q56I , p7[ p8| ] . The numerical results obtained so far show that the length scale related to magnetic 



screening is larger than the one related to electric screening. 



5. Chrial transition 



The Lagrangian of QCD has an approximate SUa(3) chiral symmetry. This symmetry is bro- 
ken in the vacuum. The chiral symmetry breaking is signaled by non-zero expectation value of 
the quark or chiral condensate, (yy) ^ in the massless limit. This symmetry is expected to be 
restored at high temperatures and the quark condensate vanishes. There is an explicit breaking of 
the chiral symmetry by the non-zero values of u,d and s quark masses. While due to the relatively 
large strange quark mass (m s ~ 100 MeV) SUa(3) may not be a very good symmetry its subgroup 



8 



Recent highlights in lattice calculations 



P. Petreczky 



SUa(2) remains a very good symmetry and is relevant for the discussion of the finite temperature 
transition in QCD. If the relevant symmetry is SUa(2) the chiral transition is expected to be second 
order for massless light (u and d) quarks belonging to the 0(4) universality class. Recent calcu- 
lations with p4 action support this pictures [^9|]. This also means that for non-zero light quark 
masses the transition must be a crossover. The later fact seems to be supported by calculations in 
Ref. [|6C|]. The Ua(1) symmetry is explicitly broken in the vacuum by the anomaly but it is expected 
to be effectively restored at high temperatures as non-perturbative vacuum fluctuations responsible 
for its breaking are suppressed at high temperatures. If the Ua(1) symmetry is restored at the same 



temperature as the SUa{2) symmetry the transition could be first order [61]. Recent calculations 



with staggered [q2|] as well as with domain wall fermions [ 23 ] suggest that Ua (1 ) symmetry gets 
restored at temperature that is significantly higher than the chiral transition temperature. 

For massless quark the chiral condensate vanishes at the critical temperature T® and is the or- 
der parameter. Therefore in the lattice studies one calculates the chiral condensate and its derivative 
with respect to the quark mass called the chiral susceptibility. For the staggered fermion formula- 
tion most commonly used in the lattice calculations at non-zero temperature these quantities can 
be written as follows: 

Xm,q( T ) = n f d %^ { qt = Xc ' Misc + Xl' co " Q = l i S i ( 5 - 2 ) 

where the subscript x = % and x = will denote the expectation value at finite and zero temperature, 
respectively. Furthermore, D q = m q ■ 1 + D is the fermion matrix in the canonical normalization and 



rif = 2 and 1 for light and strange quark. In Eq. ( |5.2[ ) we made explicit that chiral susceptibility 
is the sum of connected and disconnected Feynman diagrams. The disconnected and connected 
contributions can be written as 



,2 



x * disc = TwW t {<( Tr V) 2 > " < TrZ V> 2 } . ^ 



Xccon = -^Tr£(D ¥ 1 (x,0)D ? 1 (0,x)) , q = l,s. (5.4) 

X 

The disconnected part of the light quark susceptibility describes the fluctuations in the light quark 
condensate and is directly analogous to the fluctuations in the order parameter of an 0(N) spin 
model. The second term (Xq.con) arises from the explicit quark mass dependence of the chiral 
condensate and is the expectation value of the volume integral of the correlation function of the 
(isovector) scalar operator yy. Let me note that in the massless limit only Xl,disc diverges. 

5.1 The temperature dependence of the chiral condensate 

The chiral condensate needs a multiplicative, and also an additive renormalization if the quark 
mass is non-zero. Therefore the subtracted chiral condensate is considered 

*< {T) -(wko-2(wW (5 - 5) 
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Figure 5: The subtracted chiral condensate calculated with HISQ action in the continuum limit compared 
to the renormalized Polyakov loop (left) and light quark number fluctuation (right). The continuum results 
for A/ s and L ren have been taken from Ref. [Q, while the continuum results for %2 nave been taken from 
Ref. M. 



In Fig. ^| I show continuum results for A/ s calculated with HISQ action and compared to the 
renormalized Polyakov loop and baryon number fluctuation previously discussed in relation to the 
deconfining transition. The rapid increase in %2 happens roughly in the same temperature interval 
where Ai, s shows a rapid decrease, while it is difficult to make similar statements for L ren due to 
its very smooth behavior. However, it is clear from Fig. |] that L ren is very far from unity for 
temperatures where A/ . s is very small. 

Another way to get rid of the multiplicative and additive renormalization is to subtract the zero 
temperature condensate and multiply the difference by the strange quark mass, i.e. consider the 
following quantity 

Aq=d + 2m s r\({yy) q ^ - (W) q ,o) q = l,s. (5.6) 



The factor r\ was introduce to make the combination dimensionless. Here r\ is the scale parameter 
defined from the zero temperature static potential [Q. It is convenient to choose the normalization 
constant to be the light quark condensate for mi = multiplied by m s r\. In Fig. |6| the renormalized 
quark condensate is shown as function of the temperature for HISQ and stout actions. We see a 
crossover behavior for temperature of (150 — 160) MeV, where Af drops by 50%. The difference 
between the stout and HISQ results is a quark mass effect. Calculations for HISQ action were 
performed for m % = 160 MeV, while the stout calculations were done for the physical quark mass. 
For a direct comparison with stout results, we extrapolate the HISQ data in the light quark mass 
and also take care of the residual cutoff dependence in the HISQ data. This was done in Ref. [19] 
and the results are shown in the figure as black diamonds demonstrating a good agreement between 
HISQ and stout results. Contrary to Af the renormalized strange quark condensate Af shows only 
a gradual decrease over a wide temperature interval dropping by 50% only at significantly higher 
temperatures of about 190 MeV. The subtracted chiral condensate has also been calculated using 
Wilson and overlap formulations [2C, p4j ]. These calculations show good agreement with staggered 
results at the corresponding values of the pion mass. 
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Figure 6: The renormalized chiral condensate Af for the HISQ action with m/ / m s — 0.05 is compared to the 
stout data. In the right panel, we show the renormalized strange quark condensate Af for the HISQ action. 

5.2 O(N) scaling and the transition temperature 

In the vicinity of the chiral phase transition, the free energy density may be expressed as a sum 
of a singular and a regular parts, 



/ = ~y taZ = f s ing(t,h) + f reg (T,mi,m s ) . 



(5.7) 



Here t and h are dimensionless couplings that control deviations from criticality. They are related 
to the temperature T and the light quark mass m/ as 



1 T - T L ° 



h=—H 
ho 



H 



mi 



m, 



(5.8) 



where denotes the chiral phase transition temperature, i.e., the transition temperature at H = 0. 
The scaling variables t, h are normalized by two parameters ?o and ho, which are unique to QCD 
and similar to the low energy constants in the chiral Lagrangian. These need to be determined 
together with T c °. In the continuum limit, all three parameters are uniquely defined, but depend on 
the value of the strange quark mass. 

The singular contribution to the free energy density is a homogeneous function of the two 
variables t and h. Its invariance under scale transformations can be used to express it in terms of a 
single scaling variable 



t/h^ 8 



-T9 fho 

77 



to r« 



i/ps 



l J _ j0 / i \ i/ps 

77 



~0 



(5.9) 



where j3 and 8 are the critical exponents of the 0(N) universality class and zo = to/h^ 8 . Thus, 
the dimensionless free energy density / = / /T 4 can be written as 



f(T,m,,m s )=h l+1 / S f f (z)+f reg (T,H,m s ) , 



(5.10) 



where fj is the universal scaling function and the regular term f reg gives rise to scaling violations. 
This regular term can be expanded in a Taylor series around (t,h) = (0,0). 
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It should be noted that the reduced temperature t may depend on other couplings in the QCD 
Lagrangian which do not explicitly break chiral symmetry. In particular, it depends on light and 
strange quark chemical potentials pL q , which in leading order enter only quadratically, 




The transition temperature can be denned as peaks in susceptibilities (response functions) that 
are second derivatives of the free energy density with respect to relevant parameters. Since there 
are two relevant parameters we can define three susceptibilities: 

d 2 f d 2 f _ d 2 f 

Xm.l — ^ 9 j Xt.l — "T~=i ) Xt,t — ~T7T • (5-12) 

omf dtdmi at z 

Thus three different pseudo-critical temperatures T m j, T t j and T tJ can be defined. In the vicinity 
of the critical point the behavior of these susceptibilities is controlled by three universal scaling 
function that can be derived from ff. In the chiral limit T m j = T t j = T t t = r c °. There is, however, 
an additional complication for 0(N) universality class: while % m j and %t,l diverge at the critical 
point for m/ — > 

X m ,i~m} /5 -\ Xa- m f- m8 , (5-13) 

Xt,t is finite because a < for 0(N) models ( Xt,t ~ \t\~ a )■ Therefore, one has to consider the 
third derivative of / with respect to t : 

d 3 f 

Xw = -^- (5-14) 

In the vicinity of the critical point the derivatives with respect to t can be estimated by taking 
the derivatives with respect to jj.?, i.e. the response functions Xt.l and %t,t,t are identical to the second 
Taylor expansion coefficient of the quark condensate and the sixth order expansion coefficient to the 
pressure, respectively. The former controls the curvature of the transition temperature as function of 



the quark chemical potential \i q and was studied for p4 action using N t = 4 and 8 lattices Q63p - The 
later corresponds to the sixth order quark number fluctuation which is related to the deconfinement 
aspects of the transition. The fact that this quantity is sensitive to the chiral dynamics points to 
a relation between deconfining and chiral aspects of the transition. In the following I discuss the 
determination of the transition temperature defined as peak position of % m i, i.e. T c = T m j. 

5.3 Determination of the transition temperature 

The 0(N) scaling described in the above subsection can be used to determine the pseudo- 
critical temperature of the chiral transition. For the study of the 0(N) scaling it is convenient to 
consider the dimensionless order parameter 

M h = m s ^^. (5.15) 

The subscript "b" refers to the fact that this is a bare quantity since the additive UV divergence is 
not removed. From the point of view of the scaling analysis this divergent term is just a regular 
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contribution. For sufficiently small quark mass and in the vicinity of the transition region we can 
write 

M b (T,H)=h 1 / 5 f G (t/h l ^ s ) + f Mtreg (T,H). (5.16) 
Here fc(z) is the scaling function related to ff and was calculated for 0(2) and 0(4) spin models 



541 p5|, pq ]. The regular contribution can be parametrized as [ |19| ] 

fM, reg (T,H) = a t (T)H 



a + a l —^ + a 2 ( ) 



Then we have the following behavior for the light chiral susceptibility 

-2 



Xm,l T 2 ( 1 1/5 _! df M . reg (T,H) \ 



™& f x {z) = \\f G {z)-jf G {z)\. (5.18) 
One then performs a simultaneous fit to the lattice data for M b and Xm,l treating T® ,to,ho,ao,ai 



and (32 as fit parameters Q19Q. This gives a good description of the quark mass and temperature 
dependence of %m,i and allows to determine accurately the peak position in %m,h Using this scaling 
analysis T c has been determined for asqtad and HISQ actions for different N?. Having determined 
T c for HISQ and asqtad action for each A^ T a combined continuum extrapolation was performed 



using different assumption about the N t dependence of T c which resulted in the value [ ]19[ ] : 



T c = (159 ±9) MeV. (5.19) 

The analysis also demonstrated that HISQ and asqtad action give consistent results in the con- 
tinuum limit. The Budapest- Wuppertal collaboration found T c = 147(2)(3)MeV, 157(3) (3)MeV 
and 155(3)(3)MeV defined as peak position in the chiral susceptibility, inflection points in A/ , s and 
respectively These agree with the above value within errors. The peak position in Xdisc 
calculated using Domain Wall Fermions is also consistent with the T c value in Eq. ( |5.19 ). 



6. Conclusions 



In recent years significant progress has been made in lattice QCD calculations at non-zero tem- 
perature. Chiral and deconfining aspects of the QCD transition have been studied using improved 
staggered quark formulation allowing to control discretization effects. Some quantities have been 
calculated at small baryon density using Taylor expansion in chemical potentials. At sufficiently 
low temperatures lattice results can be understood in terms of hadron resonance gas model, while 
at high temperatures resummed perturbative calculations describe the lattice data quite well. For 
several quantities it has been shown that in the continuum limit different discretization schemes, 
including discretizations other than staggered, give consistent results. In particular, agreement has 
been reached on the value of the chiral transition temperature. There is still disagreement in the 
lattice calculation of the equation of state. 
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